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Abstract. We present the results from a multi- 
wavelength study of the mixed H+He composition DAO 
white dwarf RE J0720-318. A detailed analysis of UV 
and EUV spectroscopic data with state-of-the-art non- 
LTE photospheric models demonstrates that the observed 
opacity to EUV radiation probably results from a more 
complex structure than a simple H+He, chemically layered 
atmosphere. Instead, EUV photometry and phase resolved 
EUV spectroscopy indicate a likely spatial non-uniformity 
in the surface distribution of helium, which is consistent 
with a model in which material is accreted from the wind 
of the dM secondary. The rotational modulation of the 
spatially inhomogeneous EUV opacity allows us to esti- 
mate the rotation period of the white dwarf (0.463 ±0.004 
days). We have also reviewed two plausible origins pro- 
posed by Burleigh et al. (1997) and Dupuis et al. ( |1997a| ) 
to account for the unusual N(HI)/N(HeI)^l along this 
line of sight. We conclude that it is probably due to the 
presence of a cloud of ionized gas along this line of sight, 
rather than a circumbinary disk. The cloud, residing be- 
tween 123-170pc distant in the direction of the CMa ISM 
tunnel, may be <; 40pc in length. 

Key words: stars: abundances - stars: atmospheres - 
stars: binaries: close - stars: individual: RE J0720-318 



1. Introduction 

Soft X-ray and EUV observations of DA white dwarfs 
have revealed photospheric opacity in excess of that pro- 
vided by a pure hydrogen plasma for the majority with 
T ofr > 50 000 K (e.g. Marsh et al. |1997] ; Wolff et al |l996 ; 

1993| ) ' 



Finley 1996 ; Jordan et al. 1994 ; Barstow et al. 
In addition, high resolution IUE echelle spectra of those 
stars exhibiting the greatest levels of opacity show ab- 
sorption lines of heavy elements (e.g. C, N, O, Si, Fe 
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and Ni; Bruhweiler & Kondo 1981 ; Dupree & Ray mond 
1982t Bruhweiler & Feibelman]199l| Holberg et al. |1994 ; 
Holberg, Barstow and Sion 1998 ). In many cases these 
have been shown to be consistent with their location in 
the s tellar photosphere (e.g. G191-B2B, Reid & Wegner 
1988 ). Early theoretical work had suggested that the at- 
mospheres of DAs were more or less pure hydrogen as 
heavier elements should rapidly sink out in the strong 
gravitational field (Schatzmann 1958 ). However, further 
detailed calculations have shown that finite abundances 
of metals can be supported against gravitational settling 
by the net transfer of upward momentum from the in- 
tense radiation field to the heavy clement ions via their 



line transitions (e.g. Vauclair et al. 1979 and most recently 
Chayer et al. 1995| ). In contrast, similar calculations have 



demonstrated that helium cannot be radiatively levitated 
in observable quantities in all but the hottest objects, due 
to a relatively small number of bound-bound transitions, 
and should sink out of the photosphere on a time-scale of 



days (Vennes et al. 1988). 

White dwarfs are classified according to their opti- 
cal spectra and generally fall into one of two categories. 
The hydrogen rich DAs and the helium rich DOs dis- 
play solely strong H and He lines respectively. The DAOs 
are a hybrid class and show both strong H lines and a 
weak 4686A Hell feature. Lying in the temperature range 
50 000 - 70 000K, Fontaine & Wesemael ( |l987p proposed 
that they were transitional objects forming an evolution- 
ary link between the DOs and DAs. As a DO cooled, resid- 
ual hydrogen in its atmosphere was expected to float to 
the surface forming a layer of sufficient thickness that spec- 
troscopic traces of helium were eliminated, at least in the 
optical waveband. The white dwarf would transform first 
into a DAO, then a DA. In either case it was expected 
to have a chemically stratified atmosphere where a thin 
layer of hydrogen over-lies a predominantly helium enve- 
lope in diffusive equilibrium. At 30 000K, convective mix- 
ing would then dredge helium back into the optical photo- 
sphere, effectively returning it to the helium rich cooling 
sequence. Many observations could be accounted for with 
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this evolutionary hypothesis, e.g. the apparent absence of 
DAs with T cff > 70 000K and the DO/DB gap between 
30 000 — 45 000K, a region of the cooling sequence where 
no helium rich objects are observed. However, doubt was 
cast on this interpretation by the detailed optical study of 
a large sample of DAOs by Bergeron et al. (1994). In con- 



trast to the theoretically predicted stratified atmosphere, 
they found for many objects in their sample that the ob- 
served 4686A Hell line profile was reproduced significantly 
better by a model in which helium was homogeneously 
mixed into a predominantly hydrogen atmosphere. In ad- 
dition, the low surface gravities of many DAOs suggested 
that they were products of extreme horizontal branch evo- 
lution (EHB), inconsistent with the generally accepted 
post-AGB e voluti onary status of the DOs and DAs (e.g. 
Weidemann |l990 ). 

RE J0720-318 is one of only three DAOs detected in 
the EUV sky surveys of the ROSAT Wide Field Camera 
(WFC, Pounds et al. |l993[ ) and t he Ex treme Ultravio- 

Optical follow 



1994) 



let Explorer (EUVE, Bowyer et al 
up observations revealed a post common envelope (CE) 
DAO+dM binary with a period of 1.26 days (Vennes 
& Thorstensen Il994a| 



to speculate that the white dwarf may have been largely 
stripped of hydrogen by mass transfer during the prior 
common envelope phase. Subsequently a weak wind was 
proposed to be carrying helium up into the photosphere or 
alternatively accretion from the secondary's wind deposit- 
ing helium into the optical line forming region to produce 
the 4686A Hell feature. 

We report here the results of a new multi-wavelength 
study of RE J0720-318. Using state-of-the-art non-LTE 
models we are able to demonstrate that the opacity ob- 
served in the EUV waveband probably arises from a more 
complex structure than a simple, LTE H+He stratified 
atmosphere. Using EUV photometry and phase resolved 
spectroscopy we also show that it is likely the helium abun- 
dance is not spatially uniform across the photosphere, in 
support of a wind accretion model. As a result we have 
been able to measure the rotation rate of the white dwarf 
(0.463 ±0.004 days). 



2. Observations and data reduction. 



Barstow et al 

Thorstensen 1996). Intriguingly all three of the DAOs de- 
tected in the EUV reside in close, pre-cataclysmic variable 
(pre-CV) systems and have masses typical of isolated DAs 
(e.g. Bergeron et al. 1992j ) and consistent with post-AGB 
evolution. Radiative levitation has been shown to be in- 
sufficient to account for the observed abundance of helium 
in these objects (Vennes et al. 1988| ) and it is thought bi- 
narity is most likely responsible (Tweedy et al. 1993). It is 



1991 Vennes & 2 - L The EUVE data 



feasible that tidal effects, accretion from the secondary's 
wind and the earlier period of common envelope evolution 
may all have played or be playing a role in deter minin g 
the observed abundance patterns (Bergeron et al. 1994 ). 



Current observational evidence seems to favour a wind ac- 
cretion mechanism. For example, periodic variations in the 
soft X-ray/EUV flux fro m the DA2+KV pre-CV system 
V471 Tauri (Jensen et al. 1986 ) can be consistent with the 
rotational modulation of regions of increased helium and 
meta l opaci ty on the surface of the white dwarf (Dupuis 
et al 1997b ). Recently, time resolved studies of a Zeeman 
split 1286A Silll line believed to be formed in the photo- 
sphere of the V471 Tauri white dwarf have revealed a vari- 
ation in its equivalent width on the soft X-ray/EUV and 
optical period (Sion et al. 1998). The time of greatest line 
strength corresponds to soft X-ray/EUV minimum and 
optical maximum, more or less confirming the magnetic 
wind accretion interpretation first proposed by Jensen et 



An EUVE observation of RE J0720-318 was carried out 
during the period 1995 December 16-20 and is now avail- 
able from the public archive at Goddard Space Flight Cen- 
tre. The target was observed in dither mode, which con- 
sists of a series of discrete pointings slightly offset from 
each other and is carried out to minimize the effects on 
the data of detector spatial response variations. RE J0720- 
318 was detected in all three spectrometers and by the 
Deep Survey (DS) imager with effective exposure times 
of 122 459s (SW 70-190A), 121067s (MW 140-380A), 
112 861s (LW 280-760A) and 112 564s (DS 70-180A). 
Standard IRAF EUV and XRAY/PROS software were 
used to extract light-curves of RE J0720-318 from each of 
the four detectors, following procedures described in de- 
tail on the CEA website at Berkeley^. The count rates for 
all four were corrected for detector dead-time and an ef- 
fect known as "Primbsching" which occurs when the satel- 
lite telemetry buffer saturates, resulting in some photon 
events being lost. Additionally, the DS count rates were 
modified to account for times when the source lay in the 
deadspot, a region of the detector with reduced gain (e.g. 
Sirk et al. 1997). The spectral extractions were also per- 



al. (1986). The influence of tidal forces, the CE phase and 
weak mass-loss is less clear and observations are essential 
for constraining theoretical models of close binary evolu- 
tion. In an initial analysis of the EUVE spectrum of RE 
J0720-318, Burleigh et al. p97| ) were able to reproduce 
the data with a simple LTE H+He stratified model. The 
unexpectedly low H-layer mass of 3 • 10~ 14 Mq led them 



formed using IRAF/EUV software, following our standard 
procedures for reducing spectroscopic data which includes 
quadratically adding a 5% systematic error to the data to 
account for residual detector fixed pattern efficiency vari- 
ations (FPN). A more detailed description of these meth- 
ods m ay be found in our earlier work (e.g. Barstow et al. 



1997) 



http:/ /www.cea.berkeley.edu 
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Table 1. The HST observation log - the time and the 
date given are for mid-exposure; phase is calculated from 
the Vennes & Thorstensen ( |1996 ) ephemeris, assuming P 
= 1.26245 days and T (HJD) = 2449735.327. 



Obs. ID. 


Date 


Timc(UT) 


A(A) 


Exp.(s) 


<S> 


z3c80405 


17/09/96 


00:06 


1623-1658 


3x571 


0.74 


z3c80408 


17/09/96 


01:35 


1623-1658 


4x598 


0.79 


z3c80305 


25/10/96 


11:08 


1532-1568 


3x571 


0.21 


z3c80308 


25/10/96 


12:38 


1532-1568 


4x598 


0.26 


z3c8030b 


25/10/96 


14:15 


1623-1658 


4x598 


0.31 


z3c8030c 


25/10/96 


15:51 


1623-1658 


4x598 


0.36 



2.2. The HST GHRS and optical observations. 

A total of six HST GHRS observations of RE J0720-318 
were obtained from the public archive at STSCI. These 
were taken with the G160M grating through the Large Sci- 
ence Aperture (LSA) and cover a waveband of 35A with a 
resolution better than O.lA. Each dataset consists of sev- 
eral shorter exposures which were co-added to produce the 
resulting spectra. A Pt /Ne hollow cathode lamp exposure 
(WAVE) was acquired immediately prior to each of the 
science observations, allowing us to use IRAF/STSDAS 
routines to re-calibrate the wavelength solution and sub- 
sequently reduce the uncertainty in the wavelength scale 
from 15kms _1 to 4kms . The observations were timed to 
coincide with binary quadrature as predicted by the opti- 
cally derived ephemeris of Vennes & Thorstensen ( 199G| ). 
Details of these, including wavebands, exposure times and 
the predicted binary phase are shown in Table [l]. The anal- 
ysis of UV and EUV data is complemented by a re-analysis 
of a co-added optical spectrum of RE J0720-318 obtained 
on the 1.9m Radcliffe telescope at SAAO in February 
1995. Full d etails of this dataset may be found in Barstow 
et al. ( |1995[ ). 



3. Analysis and results. 

3.1. Models and spectral analysis techniques. 

Our non-LTE photospheric models are generated with the 
latest version of the atmosphere code TLUSTY (version 
195, Hubeny 1988; Hubeny and Lanz 1995) and the spec- 



tral synthesis program SYNSPEC (Hubeny et al. [19940 
and include the full effects of line blanketing. Detailed hy- 
drogen and helium line profiles are calculated from the ta- 
bles of Schoning & Butler (private comm.) and Scheming 
& Butler ( 1989| ) respectively. To account for the opac- 
ity to EUV radiation of HI, Hel and Hell along the line 
of sight, during the analysis of EUV spectroscopic data, 
photospheric models are used in conjunction with the ISM 
model of Rumph et al. (1994). This has been modified to 



include the effects of the converging Lyman line series of 
Hel and Hell at 504A and 228A respectively. 

Our standard spectral analysis techniques have been 
descri bed i n detail in previous publications (e.g. Barstow 
et al. 1997 ) and will only be briefly re-iterated here. We 



carry out comparisons between models and data using 
the spectral fitting program XSPEC (Shafer et al. 1991). 
XSPEC works by folding a model through the instrument 
response before comparing the result to the data by means 
of a x 2 -statistic. The best fit model representation of the 
data is found by incrementing free grid parameters in 
small steps, linearly interpolating between points in the 
grid, until the value of \ 2 is minimized. Errors are calcu- 
lated by stepping the parameter in question away from its 
optimum value and re-determining the minimum \ 2 - This 
is repeated until the difference between the two values, 
A% 2 , corresponds to la for a given number of free m odel 
parameters as prescribed by Lampton et al. ( 1976 ). All 
errors quoted in this work are la unless stated otherwise. 

3.2. The optical data analysis 

Previous studies of the optical spectrum of RE J0720-318 
have employed LTE H+He models (e.g. Barstow et al. 



1995). A re-analysis of the optical spectrum has been car- 
ried out with a grid of non-LTE, homogeneous, H+He 
models, primarily for consistency, since non-LTE models 
are used here in the analysis of both the UV and EUV 
data. Recently it has been demonstrated that estimates of 
the effective temperatures of hot DAs determined by fit- 
ting non-LTE model Balmer lines to the observed profiles 
are systematically lower by « 4000K at w 60 000K com- 
pared to the earlier LTE analyses (Barstow et al. 



1998) 



Barstow et al. ( |1998| ) als o repo rt that in contrast to the 
results of Bergeron et al. ( 1994 ) who found that the inclu- 
sion of small amounts of helium into LTE models lowers 
the determined effective tem peratu re by several 1000K at 
« 60 000K (Bergeron et al. 1994 ), no similar effect oc- 
curs in non-LTE. Fitting our grid of non-LTE homoge- 
neous H+He models we determine T cff = 56600T.^K, 
logg = 7.41 + 0.1 and a helium abundance of log (Hc/H) = 
-3.26lg;23. We adopt T cB = 55 000K and logg = 7.5 for 
the subsequent analysis of UV and EUV data. 

3.3. The GHRS spectra 

3.3.1. Line velocity measurements 

The two spectra covering the 1532- 1568 A range both dis- 
play a prominent CIV resonance doublet. A glance at 
Fig. |l| indicates that the doublet profiles may be complex 
with weaker, red-shifted absorption lines present besides 
the main components. A test of the significance of these 
was carried out by applying an F-test to the fit statistics 
obtained using XSPEC to match both single and double 
gaussian profiles to the observed lines. The result of this 
suggests they are real features (> 99.999% confidence). A 
1640A Hell absorption line is displayed in the four spec- 
tra covering the 1623-1658A waveband. This feature was 
assumed to arise in the stellar atmosphere as photospheric 
helium has previously been detected in both the opti- 
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Fig. 1. The spectrum z3c80308 covering the CIV doublet. The line profiles may be complex and display weaker 
absorption redward of the main components. 



cal (Barstow et al. 1995; Vennes & Thorstensen 1994a) The Hell velocities were used in conjunction with the 



and the EUV spectrum (Burleigh et al. 1997; Dupuis et 
al.|1997aD of RE J0720-318. Although the four spectra cov- 



ering this waveband are rather noisy, close examination of 
z3c8030b suggests that there may also be a weak redward 
component to this profile, which can be consistent in ve- 
locity space with the weaker CIV features. However, from 
fitting gaussains to similar looking features in this spec- 
trum which are attributable to noise, it was found that 
these could also be identified as absorption lines at a sim- 
ilar level of confidence (< 99%). Furthermore, the spec- 
trum taken immediately after and covering the same wave- 
band (z3c8030e), does not appear to display this feature. 
Subsequently, the velocities of the lines in the six spec- 
tra were measured using the STARLINK package DIPSO 
to fit gaussians or where appropriate, double gaussians, 
to the observed profiles. The errors were calculated by 
quadratically adding the uncertainties in both the fitting 
procedure and in the wavelength calibration. The helio- 
centrically correct velocities along with their associated 
errors are given in Table @. 



binary cphemeris of Vennes & Thorstensen (1996) to con- 
strain the origin of features observed in the 1532-1568A 
waveb and. Unfortunately, since Vennes & Thorstensen 
( 1996| ) do not give their measured 4686A Hell line ve- 
locities in tabular form, it has not been possible to fully 
update this cphemeris. Instead, the Hell measurements 
were used to adjust K^d and 7t t, the white dwarf semi- 
amplitude and the velocity offset (gravitational and sys- 
tem) respectively. Adopting the Vennes & Thorstensen 
binary period of P = 1.26245 days and epoch(HJD) of 
T = 2449735.327, both of which should be sufficiently 
accurate for the purposes here, sinusoids of the form 
V(T) = j tot + K WD sin [tt + 2tt(T - T„) / P] were fit 
to the four data points. The best match was found be- 
tween the observed and predicted velocities for a value of 
K W d = 76.2±5.9kms~ 1 and j tot = 38.4±6.0kms _1 . Com- 



paring these figures to the Vennes & Thorstensen (1996) 
values of K\vd = 114±13kms _1 and jtot — 51±10kms _ , 
within their respective errors, consistency is found be- 
tween estimates of "ftot- However, the difference in Kwd 
is rather large. It is possible that the weak 4686A Hell 
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Fig. 2. The best fit homogeneous (black) and stratified (grey) models to the GHRS spectrum z3c8030b. See Table |^ 
for details. 



profile from which their values for Kwd and "/tot were de- 
termined, may be contaminated with emission from the ir- 
radiated face of the secondary star. It i s inter esting to note 
from Fig. 2 of Vennes & Thorstensen ( 1996 ) that the two 
velocities which agree best with the lower value for Kwd 
are obtained at phases 0.783 and 0.828 where re-processed 
emission is approaching a minimum and the Hell pro- 
file appears deepest. Furthermore, Vennes & Thorstensen 
( |l996|) point out that the errors on their velocity values are 
formal and are likely to underestimate the true uncertain- 
ties. Comparing the measured CIV line velocities to the 
predictions of the modified ephemeris, it was found that 
the main components are consistent with the predicted 
photospheric velocity of the white dwarf. If this interpre- 
tation is correct then this represents the first direct de- 
tection of carbon in the atmosphere of this white dwarf. 
The weaker components, are red-shifted by ~ 50kms _1 
with respect to the apparent photospheric velocity. We 
note that the system velocity is 28 ± 3kms _1 (Vennes & 



Thorstensen 1996). The possible origin of these features 



will be discussed later. 



Table 2. The observed and laboratory wavelengths for the 
lines detected in each of the GHRS spectra. Top, the lines 
interpreted as photospheric and lower the possible weak 
red-shifted CIV components. The heliocentric velocities of 
these lines are also shown. 



Obs. ID. 


Ion 


Lab. A 


Obs. A 


Vcl. 






(A) 


(A) 


(km s _1 ) 


405 


Hell 


1640.474 


1641.099 ± 0.030 


114.2 ± 6.8 


408 


Hell 


1640.474 


1641.088 ± 0.048 


112.2 ± 9.7 


305 


CIV 


1548.195 


1548.052 ± 0.043 


-27.7 ± 9.3 


305 


CIV 


1550.770 


1550.606 ± 0.015 


-31.7 ± 5.0 


308 


CIV 


1548.195 


1548.027 ± 0.014 


-32.5 ± 4.9 


308 


CIV 


1550.770 


1550.606 ± 0.030 


-31.7 ± 7.1 


30b 


Hell 


1640.474 


1640.286 ± 0.025 


-34.6 ± 6.1 


30c 


Hell 


1640.474 


1640.379 ± 0.030 


-17.4 ± 6.9 


305 


CIV 


1548.195 


1548.282 ± 0.015 


16.8 ± 5.0 


305 


CIV 


1550.770 


1550.854 ± 0.033 


16.2 ± 4.8 


308 


CIV 


1548.195 


1548.285 ± 0.026 


17.4 ± 6.5 


308 


CIV 


1550.770 


1550.871 ± 0.030 


19.5 ± 7.6 



3.4- The abundances of He and C from UV data 

Two small grids of models were generated to measure the 
photospheric abundances of helium and carbon. The first 
grid consists of models including a homogeneous mix of 
H, He and C. The second grid consists of models with a 
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200 
Wavelength A 

Fig. 3. The non-LTE model fits to the EUV spectrum of RE J0720-318, a. the homogeneous H+He+C (dark grey), 



N(HI) = 2.5 • 10 18 cm- 2 , N(Hel) = 1.5 • 10 18 c m - 2 , N(HeII) = 7.5 • 10 17 c m - 2 , log (He/H) = -3.8, log (C/H) = -5.5, b. 
the stratified H+Hc+C (light grey), N(HI) = 2.5 • 10 18 cm- 2 , N(Hel) = 1.6 • 10 18 cm 2 , N(HeII) = 5.1 • 10 12 c m - 2 , log 
(H-laycr) = -14.15, log (C/H) = -5.6 and c. the homogeneous models incorporating additional heavy elements (black), 
see Table § 



stratified atmospheric structure where the distribution of 
helium in the predominantly hydrogen envelope was cal- 
culated by integrating the equations of ordinary diffusion 
for a given hydrogen layer mass (Vennes et al. 1988). The 



homogeneous model grid was specified by the parameters 
log (He/H) and log (C/H) and the stratified grid by the 
parameters log(H- layer mass) and log (C/H). The anal- 
ysis of the UV spectra was carried out in a three stage 
procedure which was repeated for each grid. First the car- 
bon abundance was measured allowing the helium abun- 
dance/hydrogen layer mass to settle at an arbitrary value 
between the limits of the calculated grid. Second, with the 
carbon abundance frozen at the average value from the 
first stage, the grid was fit to the four spectra displaying 
the helium absorption profile. In the final step, with the 
helium abundance/hydrogen layer mass held at the aver- 
age value obtained in the previous step, the carbon abun- 
dance was re-measured. The results from the two latter 
stages are given in Table ||. The best fitting homogeneous 
and stratified models to spectrum z3c8030b are shown in 



Fig. |^. It is clear that while the stratified structure predicts 
a Hell line core which is too shallow, the homogeneous 
model with a helium abundance of log (He/H) » —3.74 
produces a good match. Similarly, in the remaining three 
spectra covering this wavelength range, the shape of the 
Hell 1640 A line is only accurately reproduced by a homo- 
geneous model. 

3.5. The EUVE data 
3.5.1. The EUV spectrum 

The grids of models calculated for the analysis of the UV 
spectra were extended to cover the EUV wavelength range. 
During the fitting process, the synthetic spectra were nor- 
malized to a V magnitude point and the ISM model pa- 
rameters HI, Hel and Hell allowed to vary freely. The 
hydrogen layer/helium abundance and carbon abundance 
grid parameters were constrained by the most extreme la 
limits determined in the analysis of the UV data. An ad- 
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Fig. 4. Inset: The DS light-curve. Main: Power spectrum of DS light curve. Power at 14.984 days 1 and 30.007 days 1 
can be attributed to the satellite orbit. The peak at 2.160 days -1 is significant at the 3<r confidence level. 



Table 3. The best fitting homogeneous model abundances 
to the photospheric lines observed in the GHRS spectra. 
Also given, where applicable are the best fit H-layer masses 
for the stratified model structure. (* limit of calculated 
model grid). 



Obs. ID. Element 



Abundance 
log (N/H) 



H layer mass 
log (M ) 



z3c80405 


He 


1.UU_ 25 


-14.301S"m 


z3c80408 


He 


-4.02l° ;» 


-U.29t° il 


z3c80305 


c 


-5.46l° ;l 6 


not applicable 


z3c80308 


c 


r 9 o+0.12 


not applicable 


z3c8030b 


He 


o 7/1 +0.20 
J - '^-0.24 


-14.42tg;S|, 


z3c8030e 


He 


-4 04 +0 ' 22 


-i4.36±g:i* 



equate fit could not be obtained with either model. As 
shown in Fig. || the best fit stratified structure predicted 
a flux between 80-240A several orders of magnitude lower 
than observed, effectively ruling out this layered atmo- 
spheric configuration. On the other hand, although the 
homogeneous model could provide a good match to the 
observed spectrum at wavelengths above 180A, shortward 



of this the predicted flux is too high. However, after the 
inclusion of other heavy elements including N,0,Si,Ne and 
Fe in the abundances shown in Table ^, an excellent match 
between the synthetic spectrum and data was achieved. 
Given that it was necessary to include these elements, 
this analysis strongly suggests that other metals, besides 
carbon, are present in the atmosphere of RE J0720-318. 



We note that Dupuis et al. (1997a) have previously fit the 
spectrum with a homogeneous, LTE model including trace 
heavy elements. 



3.5.2. The EUV Light Curve of RE J0720-318. 

Light-curves of RE J0720-318 were extracted from the 
four detectors onboard EUVE and were analysed with a 



Fourier routine based on the Lomb algorithm (Lomb 1976) 
to look for power between 0.3 - 560.0 days -1 . The three 
spectrometer light-curve power spectra all have high noise 
levels and no significant periodicities other than those re- 
lated to the satellite orbit were found. Nevertheless, we 
noted a weak feature in the SW power spectrum centred 
on ~ 2.17 days -1 . In contrast, in the power spectrum of 
the DS light-curve we found a significant peak at a confi- 
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Fig. 5. DS and SW light curves folded on our ephemeris. 



Table 4. Parameter values for the 'best fitting' global 
model to the EUVE data. 

Model Parameter Value(±l<r ISM columns) 



T 




55000K 




log g 




7.50 




N(HI) 


2.45(2.38 


-2.51) ■ 


10 18 cm -2 


N(Hel) 


1.36(1.30 


- 1.43) ■ 


10 18 cm -2 


N(HeII) 


7.93(4.16 - 


- 11.23) 


10 17 cm- 2 


log (He/H) 




-3.77 




log (C/H) 




-5.41 




log (N/H) 




-5.82 




log (O/H) 




-5.21 




log (Ne/H) 




-6.00 




log (Si/H) 




-6.41 




log (Fe/H) 




-6.91 





dence level greater than 3cr, centered on 2.160 days -1 . We 
note that the South Atlantic Anomaly may introduce false 
peaks into the power spectra of EUVE lightcurves, princi- 



pally at 0.98 days (Halpern and Marshall 1996), but also 
including an alias at half this cycle time. Power spectra 
of the primsching and deadtime correction tables for this 
observation do indeed reveal a periodicity on 0.989 days. 
However, the weak alias we detect at 0.50 days appears to 
be distinct from our 0.463 day period. Excluding satellite 
orbital power no other notable features are present. The 
DS light-curve and relevant section of the power spectrum 
are presented in Fig. || Fitting a sinusoid to this light- 
curve, the epoch of maximum EUV count rate (70-180A) 
was found to be T (HJD) = 2450068.677± 0.006. Folding 
the four light-curves on this ephemeris it was apparent 
that while the SW flux is also modulated on this period 
(Fig. ||), within statistical errors, both the MW and LW 
fluxes are constant in time, further suggesting that the 
0.463 day period is real. Accepting this, the flux change 
is likely to be due to either variations in Hell opacity 
along the line of sight or the rotational modulation of 
spatial inhomogeneities in the photospheric helium and 
metal distribution. Two spectra corresponding to phases 
$ = 0.4 - 0.6 and $ = 0.9 - 0.1 of the folded DS light- 
curve were extracted. These are shown in Fig. ^. It can 
be seen that the variation in EUV flux actually occurs 
shortward of A ~ 240A. This observation is inconsistent 
with the line of sight explanation outlined above. In such 
a case, any change would be restricted to A < 228A. We 
suggest therefore that the flux modulation is most likely 
to be due to a 0.463±0.004 day rotation of the white dwarf 
arising from inhomogeneous surface abundances. Such in- 
homogeneities on the surfaces of the white dwarfs in pre- 



C V syste ms V471 Tauri (Jensen et al. 1986 , D upuis et 
al. |l997bj ) and RE J1016-053 (Vennes et al. |1996| ) are be- 
lieved to be the signature of accretion from the wind of 
the secondary. 



3.5.3. Phase Resolved EUV Spectra 

To constrain the nature of the surface abundance inhomo- 
geneities, the two phase resolved spectra were examined 
with another small grid of photospheric models. This grid, 
hereafter grid 1, was built around the model shown in Ta- 
ble |^ such that this model point defined the grid centre. 
Grid 1 is specified by two independent parameters, the 
helium abundance, Y and the metal abundance, Z. The 
metal abundance ratios are fixed with respect to each 
other and are identical to those of the model described 
in Table [|. Metallicity, Z, can only vary as a whole and Z 
is directly indexed to the model carbon abundance. Dur- 
ing the fitting procedure, two copies of model grid 1 were 
simultaneously fit to the two phased spectra while hold- 
ing the line of sight columns fixed at the values given in 
Table |J. In this procedure, the x 2 -statistic provides a mea- 
sure of the combined fit quality of the models to the two 
spectra. In a first step, to test the assumption that the 
spectral change was due to non-uniformities in the photo- 
spheric helium distribution the two grid metal parameters 
were coupled to each other. The best fit to the two spectra 
was found and the x 2 noted. Second, to test the assump- 
tion that the spectral change was due to non-uniformities 
in the photospheric metal distribution the grid helium pa- 
rameters were coupled to each other. Again the best fit 
to the two spectra was found and the x 2 noted. The best 
fit model parameters for the two steps are given in Ta- 
ble |], where subscripts b and d refer to spectral phases 
<j) = 0.9 — 0.1 and <f> = 0.4 — 0.6 respectively. As the two fit- 
ting steps have the same number of degrees of freedom, we 



may use the likelihood statistic of Edwards (1972), given 
by the formula Li/L 2 = exp[(xi — Xi)/2], to calculate 
that the spectral changes are ~ 100 x more likely to be 
the result of helium as opposed to metal inhomogeneities 
on the surface of the white dwarf. Indeed, no significant 
improvement in the fit is achieved by allowing the helium 
and metal parameters to vary completely independently. 
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Fig. 6. EUV spectra of RE J0720-318 corresponding to phases $ = 0.4 - 0.6 (grey) and $ = 0.9 - 0.1 (black). It can 
clearly be seen that the change in EUV flux is confined to A < 240A. This is consistent with an apparent change in 
photospheric opacity levels. 



Table 5. Best fit model parameters assuming, in turn, 
that helium and metals are homogeneously spatially dis- 
tributed across the stellar photosphere. Y is the helium 
abundance, Z is the metal abundance. For both fits v 
1 77. 



Parameter 


Value (log (N/H)) 


z b 


-5.50(-5.53 - -5.44) 


Y b = Y d 


-3.70(-3.75 - -3.65) 


z d 


-5.37(-5.42 - -5.32) 


x 2 


250.4 


Zb = z d 


-5.44(-5.50 - -5.39) 


Y b 


-3.78(-3.84 - -3.70) 


Y d 


-3.60(-3.66 - -3.54) 


x 2 


241.3 



With this in mind we have examined the possibly that 
the photospheric helium is confined to a limited region of 
the white dwarf surface. This was carried out by attempt- 
ing to fit a linear combination of grid 1 and a further grid 
made up of models representative of the radiative levita- 



Table 6. The parameters of the model representative of 
the predictions of radiative levitation theory. 



Model Parameter 


Value(±lcr) 


T 


55000 K 


log g 


7.50 


log (He/H) 


-5.00 


log (C/H) 


-5.50 


log (N/H) 


-5.50 


log (O/H) 


-5.50 


log (Ne/H) 


-6.20 


log (Si/H) 


-7.20 


log (Fe/H) 


-5.20 - -7.00 



tion calculations of Chayer et al (1995), hereafter grid 2, 
and described in Table |^. Grid 2 has only one variable 
parameter, log (Fe/H). This may vary between an abun- 
dance of —5.2, consistent with the theoretical prediction 
of Chayer et al. ( 1995| ) and a lower limit of —7.0. The 
abundances of the other elements included, He,C,N,0,Si 
and Ne, are fixed at the values given in Table | With 
the grid 1 Y fixed at the grid upper limit of —3.15, the 
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grid 1 Z fixed at —5.5, and with the grid 2 iron abun- 
dance initially fixed at the theoretical value of —5.20, by 
allowing the ISM columns and relative model proportions 
to vary freely, no acceptable fit to the integrated EUV 
spectrum could be obtained. The predicted SW flux was 
several orders of magnitude too low. Next, the fitting pro- 
cedure was repeated with a freely varying grid 2 iron abun- 
dance. An accurate match to the EUV data was found 
with log (Fe/H) = —7.0 and relative model proportions of 
~60% and ~40% for grid 1 and grid 2 respectively. The fit 
achieved here to the EUV data is as statistically significant 
as that of the model described in Table || As this linear 
combination is also as successful in reproducing the UV 
observations as our single component model, one cannot 
be preferred over the other. The EUV spectral variation 
may be reproduced by a ^20% change in the relative pro- 
portions of the models i.e. <fi — 0.4 — 0.6 - 70% grid 1 + 
30% grid 2 and <j> = 0.9 - 0.1 - 50% grid 1 + 50% grid 2. 

4. Discussion 

4--1- Is RE J0720-318 viewed through a circumbinary 
disk? 

The two previous studies of the EUV spectrum of RE 
J0720-318 both suggested that the peculiar N(HI)/N(HeI) 
~ 1 along this line of sight may be due to the presence 



of a circumbinary disk (Burleigh et al. 1997, Dupuis et 
al. 1997a), a remnant of the prior phase of common enve- 



lope evolution. In comparing the spectrum to that of the 
DA, RE J0723-277, which lies at a distance of ~ 117pc 
and is angularly separated in sky from, the ~ 187pc dis- 



tant, RE J0720-318 by 4°, Dupuis et al. (1997a) alterna- 



tively suggested that the line of sight towards the DAO 
intercepts an ionized cloud. However, neither study made 
a comparison between the columns observed towards the 
white dwarfs and the B stars eCMa and /3CMa, two other 
objects in this region of sky which have been spectroscop- 
ically observed by EUVE. Using a grid of LTE, homo- 
geneous H+He photospheric models calculated by Detlev 
Koester our measurements of the columns towards RE 
J0723-277 are N(HI)= 8.56 - 9.72 • 10 17 cm" 2 , N(HeI)= 
1.38- 1.60- 10 17 cm" 2 and N (Hell) = 0.64- 1.98- 10 | 7 cm- 2 , 
consistent with the values found by Dupuis et al. ( 1997a[ ). 
The implied ionization fractions of fjj — 0.51 — 0.76 and 
fife = 0.28 — 0.59 for H and He respectively are consistent 
with the typical values found by Barstow et al. (1997) 
in a spectroscopic EUV study of 13 DA white dwarfs. It 
has been suggested that the similarities in ionization frac- 
tions observed towards the white dwarfs in their sample 
may indicate that the majority of the absorbing material 
lies within the local interstellar cloud (LIC; Holberg et al. 
1998| ). Indeed, Holberg et al. ( |1998[ ) have found strong 
evidence supporting this from a high resolution HST 
STIS observation of the DA white dwarf RE J1032+53. 

12-10 17 cm- 2 



put at a distance of between 123-143pc. This is very sim- 
ilar to the measurement of N(HI) towards RE J0723-277 
obtained here. Although constraints on N(HcI) towards 
eCMa are not specifically given, examination of Fig. 5b 
of Cassinelli et al. ( 1995| ) suggests that the neutral he- 
lium column along this line of sight may also be consistent 
with that observed towards the DA (~ 10 17 cm~ 2 ). In fact 
high resolution observations of the line of sight towards 
eCMa obtained with the HST echelle B and G190M grat- 
ings, indicate that most of the intervening material lies 
within 3pc of t he Su n, consistent with residence in the 
LIC (Gry et al. 1995 ). Over the remaining distance, the 
mean volume density is found to be less than 2-10~ 4 cm~ 3 . 
It seems likely then, that the opacity observed towards RE 
J0723-277 also resides within the LIC. In stark contrast, 
the columns observed towards /3CMa, only 14° away in 
the sky from eCMa, but more distant at between 139- 
170pc, are N(HI)= 2.0 - 2.2 ■ 10 1 8 cm~ 2 and N(HeI)> 



1.4- 10 8 cm 2 (Cassinelli et al. 1996). This implies a ratio 



of N(HI)/N(HeI)^ 1 along this line of sight, very sim- 
ilar to that found towards RE J0720-318. A comparison 
between the column densities measured towards these two 
objects, bearing in mind N(HI)=2.38 — 2.51 ■ 10 18 cm -2 and 
N(HeI)= 1.30 - 1.43 • 10 18 c m - 2 towards RE J0720-318, 
may in fact indicate that both /3CMa and the DAO are 
viewed through the same LISM structure. The existence 
of an ionized gas cloud lying beyond the local cloud and in 
the direction of these objects is strongly supported by the 
findings of a high resolution study of the line of sight to- 
wards /3CMa (Dupin & Gry 1998). It appears that this is 



Cassinelli et al. ( [19951) report N(HI)= 7 



to- 



dominated by two ionized clouds distinct from each other 
in velocity space by lOkms -1 , and which Dupin & Gry 
refer to as C and D. Neither has a velocity consistent with 
the LIC. Using sulphur, which is only slightly depleted 
in the interstellar environment, to map the distribution 
of hydrogen gas, they calculate that the total hydrogen 
column along this line is N(H to t)= 1-7 — 2.1 • 10 19 cm~ 2 . 
This is approximately 10 x the neutral helium column ob- 
served towards RE J0720-318. If the interpretation here, 
that these objects are viewed through the same structure, 
is correct, then such a region of ionized gas, residing at 
between 123-170pc away in the direction of the CMa tun- 
nel, extends for more than 40pc in length. Further STIS 
observations of RE J0720-318 and RE J0723-277 would be 
able to confirm these conclusions. Additionally, a long du- 
ration EUVE observation of RE J0720-318 might allow a 
more accurate determination of fu e along this line of sight 
which may, in turn, shed further light on the mechanisms 
responsible for the observed degree of He ionization. For 
example, consistency with the fractions found by Barstow 
et al. (1997) may indicate that the mechanism responsible 
for the level of ijj e , operates over distances of the order 
~ lOOpc. This might pressumably argue against a model 
in which the He ionization resulted from the ambient stel- 
lar EUV radiation field impinging on these clouds. 



wards eCMa, which Hipparcos (ESA. 1997 ) measurements 
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4-2. The elemental abundances in the atmosphere of an 
accreting white dwarf 



In the Burleigh et al. Jl997| ) analysis of RE J0720-318, 
LTE H+He models were used without the constraints im- 
posed by the GHRS data. While the EUV spectrum was 
interpreted as originating from a white dwarf with a thin 
hydrogen layer of 3 • 10 _14 Mq, a product of mass transfer 
during a prior common envelope phase, it was suggested 
that the homogeneously mixed helium observed in the vis- 
ible photosphere was accreted from the red dwarf wind. 
However, the evidence from the light-curves, the phase 
resolved spectra and the multi-wavelength non-LTE mod- 
elling reported here, supports a wind accretion origin for 
all of the observed helium, from the optical through to 
the EUV. To first order, this excessive abundance of he- 
lium, <T dex above the radiative levitation prediction, is 
maintained by a balance between the accretion of helium 
into and the gravitational settling of helium out of the 
white dwarf photosphere and it is possible to obtain a 
crude estimate of the ac cretion rate by solving equation 1 
(Sion & Starrfield 1984 ). However, such an estimate may 
be rather inaccurate for several reasons. For instance, in 
our calculation it is assumed that the relative abundance 
of elements in the incoming material have solar values, 
although it is possible that a slowly rotating weakly mag- 
netic hot white dwarf preferentially accretes certain ele- 
ments (Alcock & Illarionov 198C). Furthermore, the size 
of the area of the photosphere onto which the material is 
accreted is unknown. Nevertheless, noting the following, 
it should be possible to eliminate the latter difficulty and 
obtain a rough lower limit to the accretion rate. The spec- 
troscopic signature of a given abundance of helium in the 
photosphere may be reproduced by models making vari- 
ous assumptions about the abundance and distribution of 
this helium. For example, a uniform distribution can be 
assumed over the entire area of the visible hemisphere, or 
a higher abundance can be assumed over a smaller region, 
as was shown earlier. In fact, it is found from modelling 
and equation 1 , where M a is the mass accretion rate, A, 
the accretion surface area, p the mass density at formation 
depth of the Hell line used to measure the abundance and 
v the diffusion velocity relative to the centre of mass, also 
at the line formation depth, that the minimum accretion 
rate occurs for a uniform distribution of helium. 



(He/H), _ 

(He/H) M a + Apv 



(1) 



The diffusion velocity has been calculated from equa- 
tion 1 of Vennes et al. ( 198SQ , by incorporating a diffu- 
sion coefficient determined from the tables of Paquette et 
al. ( 1986] ). Putting the following values into equation 1 , 
A« l.M0 19 cm 2 , p « 8.5-10~ 9 gcm- 3 , log (He/H)* « -3.9 
and v w Orchis" 1 a lower limit of M a > 1 ■ 1O~ 18 M 
yr -1 is determined. In fact optical data taken around 1995 
January-February suggests that RE J0720-318 undergoes 



periods of more intense accretion (Finley et al. 1997). A 
H-layer mass of only 10 _13 Mq is capable of hiding traces 
of the underlying helium envelope in a 60000K hydrogen 



rich white dwarf (Vennes & Fontaine 1992), and it is found 
that this DAO may have acquired ~ 1.3 • 10 _12 Mq of hy- 
drogen since it emerged from the CE some 2.6 • 10 6 years 
ago. We conclude, that although RE J0720-318 most likely 
has a hydrogen layer mass greater than 10 _13 Mq now, it 
may have been almost completely stripped of a hydrogen 
layer when it emerged from the CE phase of evolution as 
initially proposed by Burleigh et al. ( 1997 ). 

Both successful model fits to the EUV spectrum of RE 
J0720-318 achieved here require a global iron abundance 
of log (Fe/H)~ —7.0, although due to the limited wave- 
length coverage of the high resolution UV data it is not 
possible to exclude the likelihood that larger amounts of 
iron are in fact present in the atmosphere, and the abun- 
dance of the other elements such as O and N is lower. 
However, by assuming that an edge or line blend would 
be visible if predicted to be larger than twice the S/N of 
the data, the EUV spectrum can be used to set an up- 
per limit of log (Fe/H) —6.2 on the iron abundance. In 
the atmosphere of an accreting white dwarf such as RE 
J0720-318, where the elemental diffusion time-scales are 
very much shorter than the evolutionary time-scale, the 
elemental abundances should reach an equilibrium some- 
where between those of the incoming material and the 
radiative levitation values (e.g Vennes et al. 1996). The 



theoretical radiative levitation prediction (Chayer et al. 



1995) which is highly sensitive to the effective tempera- 



ture and surface gravity of a star, for the abundances of 
iron and carbon is log (C/H) = —5.5 ±0.1 and log (Fe/H) 
= — 5.1 ± 0.1. Again assuming the abundance of elements 
in the accreting material to be roughly solar it would be 
expected that the observed global iron abundance should 
be at least log (Fe/H) ~ —5.2. In contrast to the observed 
log (C/H), the measured iron abundance for RE J0720- 
318 does not conform to a simple equilibrium prescription. 
Such inconsistencies between the theoretical and observed 
abundances of elements have been widely reported in pre- 
vious studies of hot white dwarfs. For example Barstow 



et al. (1996) find a silicon abundance in GD394 at least 
a factor 10 above that expected for a 38 500K DA white 
dwarf. In contrast, the silicon abundance of log (Si/H) 
= -7.98 ± 0.07 in RE J1614-085, a DA of similar T cS , is 
around a factor 10 below the theoretical value, while ni- 
trogen is found to be over abundant by a factor 1000 (Hol- 
berg et al. |1997| ). Holberg et al. Q1997Q suggest that both 
accretion and mass-loss are perhaps responsible for the 
non-equilibrium abundances observed in RE J1614-085. 
A mass-loss of only 10 _16 Mq yr _1 has been shown to be 
capable of depleting the photosphere of a white dwarf of 
silicon on a time-scale of several thousand years (Chayer 
et al. p.997|). In fact a star of RE J0720-318's metallicity 



(l/100th solar) is expected to unde rgo mass-loss at a rate 



of - 10~ 14 M 



y r " 



(Abbott |1982|) . However, any mass- 
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loss would also reduce the efficiency of accretion (Mullan 
et al. 1992). The Bondi-Hoyle formalism can be used to 
provide a theoretical estimate of the rate of spherically 
symmetric accretion of wind material. The prediction is 
extremely sensitive to the relative velocity of the incoming 
material, being inversely proportional to its fourth power. 
If the separation between the primary and secondary bi- 
nary components is taken to be ~ 4- 10 11 cm, based on the 
ephemeris of Vennes & Thorstensen ( |1996| ) , and conserva- 
tive estimates for the red dwarf mass- loss rate and relative 
wind velocity of 3 ■10~ 1 4 M ( T) yr _1 and ~ 550km s _1 re- 



spectively (Mullan |1996|) are adopted, an accretion rate of 



~ 1 • 1CP 16 Mq yr _1 is predicted. Although the size of the 
area of the photosphere accreting is not known, detailed 
non-LTE modelling of the 1640A line profile with linear 
combinations of pure hydrogen and mixed H+Hc composi- 
tion synthetic spectra, allowing the helium abundance and 
relative model proportions to vary freely, suggests that it 
is sizable. In fact if the helium abundance is increased 
above log (He/H)> —2.1 corresponding to ^ 35% of the 
visible surface, then a deterioration in the goodness of fit 
occurs at the 3<r level. A solution of equation 1 for these 
parameters has M a ^ 3 • 10~ 17 MQyr -1 . We stress that 
this estimate is extremely crude. However, in view of the 
observed upper limit on the iron abundance we suggest 
that, while material is being accreted into the atmosphere 
of this white dwarf, a weak radiatively driven wind may 
selectively expel heavy elements from the photosphere. 

A EUV study of another DAO+dM system, RE J1016 - 
053, detected in the WFC survey (Pounds et al. |l993|) 
an d hav ing a binary period of 0.789 days (Thorstensen et 
al. 1996 ) also found low globally averaged heavy element 
abundances on the white dwarf. The line of sight neutral 
hydrogen column towards this binary of 2.65 • 10 19 cm 



(Vennes et al. 1996| ) severely attentuates the EUV spec- 
trum longward of 180A making it more or less impossi- 
ble to differentiate between helium or metal spatial abun- 
dance inhomogeneities as responsible for the 30% peak to 
peak variation on a 57.3 minute period observed in the 
EUV flux. The low global heavy element abundances led 



Vennes et al. ( p_996| ) to suggest that perhaps metals were 
confined to a small region of the photosphere and were re- 
sponsible for this variation. However, modeling indicates 
that the spectral changes observed in RE J1016-053 can be 
reproduced by an apparent global change in helium abun- 
dance of ~ 0.15 dex, similar to the result obtained here for 
RE J0720-318. It seems more plausible that the observed 
spectral changes are due mainly to inhomogenieties in the 
spatial distribution of photospheric helium. 

4-3. Where do the weak red-shifted components originate? 

Weak, non-photospheric, high ionization features such as 
the resonance lines of CIV, SilV and NV have been re- 
ported in the UV s pectra of seve ral w hite dwarfs (e.g. 
Shipman et al. 1995 ; Holberg et al. 1995 ). These are often 



observed to be blue shifted with respect to the photo- 
spheric velocity and hence are frequently attributed to a 
shell of photoionized circumstellar gas possibly formed as 
a result of weak mass-loss. However, the weak CIV lines 
detected in the GHRS spectra of RE J0720-318 are red- 
shifted by ~ 50km s _1 with respect to the observed photo- 
spheric velocity, so this material must be moving towards 
the white dwarf, arguing against a circumstellar origin. 

The close proximity of e CMa and (3 CMa, separated 
in the sky from RE J0720-318 by - 6° and - 20° respec- 
tively allows an examination of the possibility that these 
features reside in the LISM. Distances to both stars are 



accurately known from Hipparcos (ESA. 1997 ) observa- 
tions, e CMa lies between 123-143pc away and j3 CMa be- 
tween 139-170pc. In their hi gh re solution study of e CMa 
with the HST, Cry et al. ( |l995| ) tentatively detect two 
separate non-photospheric components to the CIV reso- 
nance line profiles. The first, with a column density of 
4.2± 1.0- 10 12 cm~ 2 and velocity of-10±2kms~ 1 they as- 
sociate with a warm cloud probably lying beyond the LIC. 
The other, with a column desnity of 3.2 ± 1.0 • 10 12 cm -2 
has a velocity of 17 ± 2kms~ 1 which is consistent with 
the projected velocity of the LIC along this line of sight. 
The velocity and strength of the possible non-photospheric 
CIV resonance lines observed in spectrum of (3 CMa are 
also consistent with a formation in the LIC (Dupin & Cry 



1998). However, for the CIV features in the GHRS spectra 



of RE J0720-318, although the weighted mean velocity of 
17.li2.9kms" 1 suggests a formation in the LIC, assuming 
the lines to lie on the linear part of the curve-of-growth, a 
lower limit to the column density of N(CIV)<; 1 ■ 10 13 cm -2 
can be set. This may be rather higher than expected from 
the e and /3CMa observational results and the predictions 
of theoretical studies of the conductive interface where this 
species is believed to re side ( N(CIV)« 3.1±1.0-10 n cm" 2 ; 
Cheng and Bruhweilcr 1990[ ). 

Alternatively, it may be that they are related to the 
binary itself. Indeed, the radiation field from a hot hy- 
drogen rich white dwarf is capable of creating and main- 
taining highly ionized species in a gas in close enough 
proximity (Dupree & Raymond 1982). For example, non- 
photospheric CIV absorption lines have been observed at 
a wide range of binary phases in IUE echelle spectra of 
the DA+dM pre-CV binary Feige 24 (Dupree & Raymond 



1982). Vennes & Thorstensen (1994b) used a simple pho 



toionization model in their study of low density gas in 
the environment of the DA, and found for plausible elec- 
tron temperatures in the range 10 4 — 10 5 K, that while 
the majority of the total carbon would be in the form of 
CV, a fraction of between 10~ 2 — 10 -4 would be in the 
form of CIV in its groundstate. They concluded that a 
likely origin for the absorption features was in circum- 
stellar material possibly a remnant left over from the CE 



phase (Vennes and Thorstensen 1994b). As the HST ob- 



servations of RE J0720-318 were aquired close to binary 
quadrature, the possibility that they form in a cloud of 
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remnant CE gas located at the LI Lagrangian point can be 
discounted. Nevertheless, using the photoionization code 



XSTAR (Kallman & Krolik fL997|) it is found that for gas 
with a plausible 10 4 K, the elements of which are in 
solar proportions, lying outside the binary system at a dis- 
tance of ~ 1000 x the binary separation, where effectively 
all CIV is in groundstate, the fraction of the total carbon 
in the CIV state, can be between 1 • 10~ 4 — 0.5, for hydro- 
gen densities in the range 0.1 — 10 4 cm~ 3 . In contrast, the 
fraction of helium in the neutral state is 1 • 10~ 12 , with 
the majority being twice ionized. Thus it seems plausible 
that our CIV features may also originate in material ly- 
ing outside of the binary system, possibly remnants of the 
CE. Indeed, as poi nted o ut by both Burleigh et al. ( 1997 ) 
and Dupuis et al. ( 1997a ), a circumbinary disk is theoret- 
ically predicted to form from the remains of the common 
envelope gas. 

Further modelling of gas in the vicinity of the DAO 
was carried out to investigate whether these features may 
alternatively reside in a wind from the red dwarf com- 
panion. Departures from groundstate are likely negligi- 
ble at the den sities and temperature considered (Kall- 
man & Krolik 1997 ), except within a few stellar radii 
of the white dwarf where photoexcitation may be signif- 
icant. XSTAR was run to calculate the fraction of total 
gas in the form of CIV for a range of mass-loss rates 
(10~ 9 — 10 _15 MQyr _1 ). The wind density was assumed 
to drop off as the square of the distance from the source 
and the wind velocity taken as ~ lOOkrns -1 , roughly the 
difference between the radial velocity of the CIV lines and 
the ephemeris prediction for the red-dwarf velocity. At a 
series of logarithmically spaced distances from the white 
dwarf, the local number density of CIV was calculated. By 
integrating along the line of sight an estimate for the CIV 
column was obtained. For the chosen mass-loss rates the 
predicted column was found to lie in the range N(CIV)~ 
1 ■ 10 17 - 3 • 10 7 cm~ 3 . A mass-loss rate of lO^ 11 M yr- 1 
should produce N(CIV)~ 1.7 • 10 13 cm -3 consistent with 



the observational lower limit. Mullan (1996) reports that 
a detection of the inactive Barnard's Star at A = 6cm 
has been interpreted as free-free emission originating in a 
mass-loss wind of 2 • 10~ 11 MQyr _1 . Hence, it also seems 
plausible that these absorption features may be formed 
in a photoionizcd wind from the red dwarf. Furthermore, 
at the observed velocity, wind material would be moving 
towards the white dwarf and possibly accreting onto it, 
consistent with the wind accretion model cited to explain 
the EUVE data. 



5. Summary 

It has been demonstrated that all of the helium in the 
optical, UV and EUV spectra of RE J0720-318, and most 
probably heavy elements too, are likely accreted into the 
atmosphere from the wind of the dM companion. This 
conclusion is supported by the non-uniform spatial distri- 



bution of EUV opacity in the white dwarf photosphere. It 
has been found that even at the lower limiting rate of ac- 
cretion, the white dwarf may well have gained its entire ob- 
servable atmosphere from the wind of the secondary in the 
time since it emerged from the CE phase of evolution. Fur- 
thermore, the non-uniformity has allowed an estimate of 
the rotation period of the white dwarf (0.463 ±0.004days) 
to be made. The photospheric iron abundance is at least a 
factor 10 lower than predicted by a simple balance between 
gravitational settling and the opposing upward pressure of 
the radiation field. Along with the rather lower than ex- 
pected accretion rate, this has been interpreted as due to 
the possible presence of a selective mass-loss wind. Sev- 
eral possible origins for weak non-photospheric CIV lines 
observed in GHRS spectra of RE J0720-318 have been ex- 
amined. Although the line velocities are consistent with a 
formation in the LISM, the column density may be rather 
larger than expected. Alternatively, we have shown that 
this CIV may reside in photoionized gas in close proximity 
to the white dwarf. This may take the form of either resid- 
ual CE material surrounding the binary or in a putative 
wind from the secondary star. Finally, circumstantial evi- 
dence has been presented which argues against a circumbi- 
nary disk as responsible for the unusual N(HI) /N(HcI)~ 1 
along this line of sight. Instead, it has been shown that it 
is more likely due to an extended region of ionized gas 
lying in the direction of the CMa ISM tunnel. 
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